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Abstract. We observed the upper hne structure line of neutral carbon, Cl( 3 P2 — >3 Pi) (^Vest = 809 GHz), the 
12 CO(J=3— >2) line (v> le st = 345 GHz) and the 1.2 mm continuum emission from H1413+117 (Cloverleaf quasar, 
z = 2.5) using the IRAM interferometer. Together with the detection of the lower hne structure line (Barvainis 
et al. 1997), the Cloverleaf quasar is now only the second extragalactic system, besides M82, where both carbon 
lines have convincingly been detected. Our analysis shows that the carbon lines are optically thin and have 
an excitation temperature of T CX «30K. CO is subthermally excited and the observed line luminosity ratios are 
consistent with n(H2) ~ 10 3-4 cm -3 at T^ n — 30-50 K. Using three independent methods (C I, dust, CO) we derive 
a total molecular gas mass (corrected for magnification) of M(H/2) ~ 1.2 ± 0.3 x 10 10 Mq. Our observations suggest 
that the molecular disk extends beyond the region seen in CO(7-6) to a zone of more moderately excited molecular 
gas that dominates the global emission in Cl and the low J CO lines. 

Key words, galaxies: formation - galaxies: starburst - galaxies: high-redshift quasars: emission lines - quasars: 
individual H1413+117 - cosmology: observations 



1. Introduction 

Detections of large amounts of dust and gas in distant 
quasars have opened the possibility to study molecular 
gas properties in the early epoch of galaxy formation and 
to provide fundamental constraints on galaxy evolution. 
To date, CO emission has been detected in more than a 
dozen quasar host galaxies with z > 2, the record holder 
being at z = 6.4 (Walter etal. [2003, Bertoldi etal. E003)- 
Molecular gas masses in excess of 10 10 M have led to the 
hypothesis that the tremendous far-infrared (FIR) lumi- 
nosities (> 10 12 Lq) of these objects are not only powered 
by active galactic nuclei (AGN) but also by major star- 
bursts which might be forming cores of elliptical galaxies 
or bulges of massive spiral galaxies. 

Due to its strong magnification by gravitational lensing 
H1413+117 is a relatively strong emitter in various molec- 
ular lines and therefore a prime target to study the physi- 
cal conditions of the molecular gas at high redshift in great 
detail. It has been detected in the CO 3-2 (Barvainis et al. 
IP3M1 Wilner et al. HMTjI Alloin et al. ITW71 Barvainis et al. 
ITW71 B97 hereafter), 4-3, 5-4 (B97) and 7-6 (Alloin et al. 
1997J Yun et al. I1997|) transitions, the lower fine structure 



line of atomic carbon, C i( Pi 
(B97, Solomon et al. l2UU5|l . 



3 P ), (B97) and in HCN 



We use H = 75 km s 1 Mpc 1 and q 



= 0.5 



In this letter we report on the detection of the upper fine 
structure line of neutral carbon, Ci( 3 P 2 ^ 3 Pi), which adds 
an important piece of information: unlike CO, cool C i can 
be described by a 3-level system. This allows us to derive 
its excitation and to constrain the physical gas conditions 
by observations of the Ci( 3 Pi -> 3 P ) and Ci( 3 P 2 -^ 3 Pi) 
transitions only. We compare our results from C i with 
estimates based on thermal dust emission and CO. 

2. Observations and Results 

Observations were carried out with the IRAM Plateau 
de Bure interferometer during 2 nights in March 2003. 
The dual frequency setup was used to observe the 
12 CO(J=3^2) and Ci( 3 P 2 -^ 3 Pi) transitions towards 
H1413+117 (RA 14:15:46.23, Dec 11:29:44.0; J2000) si- 
multaneously. The receivers were tuned in single sideband 
mode (SSB) at 97.191 GHz and in dual sideband mode 
at 224.478 GHz/227.478 GHz. We used the standard D 
configuration (6 antennas) which results in synthesized 
beams of 7'.'5 x 4'/4 and 3'.'5 x l'/7 for the 3 mm and 1 mm 
bands respectively. At these resolutions the source is un- 
resolved. Typical SSB system temperatures were «120K 
and ~ 300 K at 3 mm and 1 mm. Flux and passband cal- 
ibration were obtained on MWC349. The nearby sources 
1413+135 and 1423+146 were used as secondary flux and 
phase calibrators. Wc estimate the flux density scale to be 
accurate to about 10%. 
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Fig. 1. Spectra of the 12 CO(J=3^2) and Ci( 3 P 2 -^ 3 Pi) 
lines towards H1413+117, superposed on their Gaussian 
fit profiles (see Table ^ for parameters). The velocity 
scale is relative to the tuning frequencies of 97.191 and 
227.478 GHz corresponding to z = 2.5579. 
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Fig. 2. Integrated Ci( 3 P 2 
in H1413+117. Contours 
and from 1.0 to 4.0 in steps of 0.5 Jybeam -1 kms" 1 . The 
rms noise is 0.5 Jybeam -1 kms -1 , the peak flux density 
is 4.3 Jybeam -1 kms -1 . The continuum has been sub- 
tracted. 



The data were processed to give data cubes with a velocity 
resolution of 10 kms -1 (3.24 MHz) at 3mm and 40 kms -1 
(29.7 MHz) at 1mm. A 3 mm continuum image was gen- 
erated by averaging the line emission-free channels. At 
1 mm the continuum image was computed from the image 
sideband. From the 1 mm signal sideband data we gener- 
ated a continuum free Ci( 3 P 2 -^ 3 Pi) spectral line cube by 
subtracting the 1 mm continuum image. The final results 
are presented in Figs.Q]&|21 

The 12 CO( J=3->2) transition {v rcst = 345.796 GHz) is de- 
tected with high signal to noise (Fig.^). Line parameters 
for the 12 CO(J=3^2) transition are similar to those re- 
ported by Wilner et al. I|1995[) and B97. The high quality 
of our 12 CO(J=3^2) spectrum allows us to determine 
the redshift of the molecular gas (as traced by CO) with 
high accuracy: z = 2.55784 ± 0.00003. 
The Ci( 3 P 2 -» 3 Pi) line (i/ rest = 809.342 GHz, Muller etal. 
200 ljl is clearly detected with a peak flux density of 
S v = 13.2 ± 2.9 mJy (Figs.^feEJ). Due to the limited band- 
width of 512 MHz (670 kms -1 ) our observations lack the 
red line-wing. However, this does not affect the determina- 
tion of the peak line intensity since the zero power level is 
well defined in interferometric observations. Tabled sum- 
marizes the line parameters derived from Gaussian fits. 
The Ci( 3 P 2 -» 3 Pi) line (z = 2.5585 ± 0.0001) is displaced 
by ps+60 kms -1 relative to 12 CO(J=3— >2) - opposite to 
the displacement of the Ci( 3 Pi— > 3 Po) transitions found by 
B97 (w-42 kms -1 ). While the difference between the CO 
and C i( 3 Pi ^ 3 Po) redshift might be attributed to low level 
baseline instabilities in the 30 m data and the poor signal 
to noise ratio, the reason for the shift in the Ci( 3 P 2 ^ 3 Pi) 
line is unclear. Gravitational amplification should not al- 



ter the center frequency for C i unless its distribution is 
different from CO, which we consider unlikely (see below). 
The apparent line shift needs to be confirmed by higher- 
sensitivity observations. If confirmed, this shift is most 
likely due to an opacity effect. 

For the dust flux, we derive an upper limit for the 
3 mm (A rcs t = 870 /im) continuum of S97.2GHZ < 1-5 mJy 
(3 c). For the 1mm (A res t = 375 /im) continuum we find 
5 , 224.5GHz = 7.5 ± 0.6 mJy. We obtained an additional 
250 GHz flux measurement using the MAMBO 117 chan- 
nel bolometer array at the IRAM 30 m telescope with 
•5250 GHz = 16 ± 2 mJy, consistent with a previous measure- 
ment (Barvainis et al. I1995f) . The radio-IR spectral energy 
distribution (SED) of H1413+117 is shown in Fig.|3 

3. Analysis 

3.1. Atomic Carbon 

3.1.1. Opacities and excitation 

Recent studies of the Milky Way and nearby galaxies 
showed that C 1 is closely associated with the CO emission 
independent of the environment (see e.g. Bennett et al. 
IT59H Ojha et al. USUI for the MW and Gerin & Phillips 
2000J Israel & Baas 120021 for nearby galaxies) . Since the 
critical density for the Ci lines and 12 CO(J = l— >0) are 
both n CI ps 10 3 cm -3 this finding suggests that the transi- 
tions arise from the same volume and share similar ex- 
citation temperatures (e.g. Ikeda etal. I2002f) . This al- 
lows us to estimate the opacities of the C 1 transitions 
via a LTE analysis assuming optically thick emission 
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Table 1. Observed line parameters. 



Line 


fobs 




AVfwhm 


I 


V s - 


L'/10 10 


Ref. 




[GHz] 


[mJy] 


[kms' 1 ] 


[Jy kms -1 ] 


[kms^ 1 ] 


[K kms" 1 pc 2 ] 




12 CO(J=3^2) 


97.1928 


30.0 ± 1.7 


416 ± 6 


13.2 ± 0.2 


-5.6 ± 2.4 


18.0±0.3 


this paper 


Cl( 3 P 2 ->- 3 Pi) 


227.4376 


13.2 ± 2.9 


368 ± 25 


5.2 ± 0.3 


53 ± 20 


1.30±0.07 


this paper 


Ci( 3 Pi^ 3 P ) 


138.351 


7.7 ± 0.8 


430 ± 46 


3.6 ± 0.4 


-47 ± 24 


2.4±0.3 


B97 



Quoted errors are statistical errors from Gaussian fits. Systematic calibration uncertainty is 10% 
a The velocity is given relative to z=2.5579, the mean redshift of the CO lines determined by B97. 

CO fluxes for 3-2, 4-3, 5-4, 7-6 are: 13.2, 21.1, 24.0, 47.3 [Jy kms" 1 ]; all values except for 3-2 are from B97. 
CO line luminosities L' are: 18.0, 16.1, 11.7, 10.9 [10 10 L Q ] 



in 12 CO(J=1^0). Since the 12 CO(J=1^0) transition 
has not been measured in the Cloverleaf we here adopt 
-k'i2co(J=i-^o) = i- 1 ^' 12 co(./^3-+2) ( see our discussion on 
the CO lines). Using the values for the line luminosi- 
ties (L') as listed in Table ^ an d assuming as a start 
T ex — Tdust = 50 K (see below) we find that both carbon 
lines are optically thin (t w = 0.14, r 2 i =0.10). For a mod- 
erate CO opacity (as suggested by B97) optical depth 
for the carbon lines are even lower. Low optical depths 
for the C i transitions are in line with findings in nearby 
galaxies (e.g. Israel & Baas 2002J ■ As discussed in Stutzki 
etal. I|1997[) the integrated line intensity ratio i?ci be- 
tween the Ci( 3 P 2 ^ 3 Pi) and Ci( 3 Pi -^ 3 Po) transition can 
be used in the optically thin limit to determine directly the 
Ci excitation temperature via T G x = 38.8 K/ln[2.11/i?ci]- 
This equation is valid if the levels are thermally pop- 
ulated. With i? C i = 0.54 ±0.07 (see Table HJ) we find 
T ex = 29 ± 3 K. For this T ex the C i lines are still optically 
thin (no = 0.15, r 2 i=0.11). 

3.1.2. Neutral carbon and molecular gas mass 

The relation between the integrated Ci( 3 P2^ 3 Pi) bright- 
ness temperature and the beam averaged C i column den- 
sity in the optically thin limit for local galaxies (z = 0) is 
given by 

where Q(T GX ) = 1 + 3e~ Tl / T - + 5e~ T2 / T - is the Ci par- 
tition function. Ti=23.6K and T2 = 62.5K are the en- 
ergies above the ground state. Correcting for redshift 
(j T m bdu — > (1 + z)~ l J T m b dv) and accounting for the 
area of the emission region (f2 s »b D\) yields the total 
number of carbon atoms in the source. Here f2 s *b denotes 
the solid angle of the source convolved with the telescope 
beam and -Da is the angular size distance to the source. 
Using the definition of the line luminosity (e.g. Solomon 
et al. I1997J) we can derive the C i mass via 

Mci = C mci Q(T CX ) I e^ T ~ L' cl{3p ^ 3Pi) (2) 



where mci is the mass of a single carbon atom and C is 
the conversion between pc 2 and cm 2 . Inserting numbers 
to Eq.H yields 

Mci = 4.566 x HT 3 Q(T CX ) i e 62 5 / T - L CI(3 p 2 ^3 Pl) [M ](3) 

For T cx = 30 K the mass of neutral carbon thus amounts 
to Mci = 2.6 x 10 7 m _1 M where m is the magnification 
by the gravitational lens. 

Deriving the mass of molecular hydrogen requires an es- 
timate of the carbon abundance relative to H2. We here 
use M 82 as a template to estimate this number. Using 
the results from Stutzki etal. ljTrj97jl . Walter etal. ©02) 
and Weifi etal. pUtjT|) we obtain [C i]/[H 2 ] « 3 x 10" 5 . 
Applying this number to the carbon mass of H1413+117 
results in M(H 2 ) » 1.4 X lO 11 m- 1 M . 

3.2. Dust emission 

The mm-wave continuum measurements combined with 
sub-mm, ISO and IRAS data (for references see cap- 
tion of Fig.BJ) are well fitted with a 2-componcnt dust 
model. To account for the steep increase of the observed 
flux with frequency in the observed mm regime we used 
K d ( t / 1 .) = 10(f r /lTHz) /3 cm 2 g- 1 with /3 = 2 (Downes & 
Solomon I1998f) . For the source solid angle of the lensed 
dust distribution we used f2g = ir m 9 2 / '4 with an angu- 
lar size of = 0.17" and a magnification of m = ll for 
both components (Venturini & Solomon E003). Masses 
and dust temperatures for both components were de- 
termined by fitting those 4 free parameters to the 12 
data points above the 97.2 GHz upper limit in Fig.|3 
For the cold dust component we find T co id rj 50 (± 2) K 
and M co id « 6.7 x 10 8 mT x M (± 12%). The warm com- 
ponent is characterized by T warm w 115 (± 10) K and 
M warm w3.9 x 10 6 m- 1 M (± 60%). This implies L Fm » 
2.4 x 10 13 mT 1 L (Lfir as defined by Helou et al. ITM^ . 
The fit results are displayed in Fig.|3 Parameters de- 
rived for the cold gas component are in agreement with 
previous studies (e.g. Hughes etal. ^97, Bcnford etal. 
1999). Using a gas to dust mass ratio of 150 (e.g. Dunne 
et al. I2000fl the total H 2 mass in the cold component is 
M(H 2 ) ~ 1.0 x 10 11 mT 1 M in agreement with the mass 
determined from Ci. Note that a larger angular size of 
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Fig. 3. Radio-IR SED of H1413+117. Displayed flux den- 
sities were taken from Barvainis & Lonsdale 1997 (20, 6, 3, 
2cm, stars), Solomon etal. 120031 (1.2 cm cross), this work 
(3mm, 1.3mm, 1.2mm, filled squares), Barvainis etal. 

(l-2mm, 760^m, 438^im, 345^m, 100/xm, 60^m, 
circles), Barvainis & Ivison 2002 (850 fim, filled trian- 
gle), Hughes etal. 119971 (800 /xm, filled circle), Benford 
et al. 119991 (350 /im, pentagon) and Rowan-Robinson 2000 
(200 /urn, triangle). The dashed lines show the thermal dust 
continuum emission for a 50 and 115 K dust component. 
The solid line is the total emission from both components. 



the cold component leads to lower temperatures for T co id, 
similarly a smaller size for the warm component will yield 
higher T warm . 

3.3. CO lines 

Our 12 CO(J=3^2) flux density is about 30% higher 
than the value reported by B97. The flux difference has 
important implications for the interpretation of the CO 
V ratios: now the 3-2, 4-3, 5-4, 7-6 CO line luminosities 
(see footnote in Tabled are decreasing with increasing J. 
So the data suggest CO is subthermally excited, at least 
for J > 3. The minimum kinetic temperature therefore 
can be lower than the lower limit of 60 K given in B97. 
We have reanalyzed the CO ratios using our 3-2 flux 
and taking the background radiation of Xcmb=9.6K 
into account. Due to the high 7-6 line luminosity relative 
to 5-4 not all CO lines can be fitted with a single gas 
component LVG model. If only the 3-2, 4-3 and 5-4 lines 
are taken into account LVG solutions can be found for 
2\in > 10 K. Since the kinetic temperature is expected to 
be close to the excitation temperature of neutral carbon 
(Israel & Baas l2TJU2l we fixed the T kin to 30 K. In a LVG 
space of [CO]/grad(V) from 1 x 10~ 5 to 2 x 10~ 4 and 
log(riH 2 ) 1.8 to 5.5, the allowed H 2 density ranges from 



3.2 < log(nH 2 ) < 4.0 with the higher limit corresponding 
to the lowest abundance per velocity gradient and vice 
versa. Independently of the abundance per velocity 
gradient or the density we find that for Tun = 30 K 
the 12 CO(J = l— +0) transition will be brighter than 
i2 CO (J=3^2) with L' 12CO(M) /L[ 2CO{J=1 ^0,9, 
consistent with the lower limit of 0.79 determined 
by Tsuboi et al. ©99 ) . Using this line ratio and 
ilW£co =0 - 8 M O (Kkms^pc 2 )' 1 (Downes & 
Solomon RWty we get M(H 2 ) » 1.6 x 10 u mT x M 
in agreement with the result for Ci and the dust. Using 
Tkin = ?dust = 50 K does not alter the results significantly. 
The predicted 7-6/3-2 ratio ranges between 0.05-0.1, 
much lower than the observed ratio of 0.6. A large 
fraction of the 7-6 line may therefore arise from much 
warmer regions such as the warm dust component. This 
is also supported by a minimum CO(7-6) excitation 
temperature of « 75 K that can be derived from the 
observed 7-6 brightness temperature. We do not think 
that selective magnification would significantly alter this 
analysis, because detailed models of the Cloverleaf lens 
show that the lens magnification is insensitive to the 
size of the mm-line emitting region (Venturini, private 
communication) . 

4. Discussion 

The mass estimates from Ci, dust, and CO give remark- 
ably similar results given that each method uses indepen- 
dent assumptions (carbon abundance, gas to dust mass 
ratio, CO conversion factor). Correcting for magnification 
we find M(H 2 )«1.2±0.3X 10 10 M Q . 
Compared to the physical conditions of the molecular 
gas determined in other dusty quasars (Tki n ~ 50-70 K, 
n(H 2 )«10 4 , see e.g. Downes etal. [1999, Carilli etal. 
2002) and local ULIRGs (Downes & Solomon H358T) the 
values determined above for H1413+117 argue for a more 
diffuse and cooler gas component. Judging from the ex- 
citation temperature of C i, the excitation of the molec- 
ular gas is between conditions found in M82's center 
(r ex ,ci= 50 K) and the mean value over the Galactic 
Center (T OX! ci = 22K) (Stutzki etal. ITMTjl . The contin- 
uum to line ratio of Lfir/£'co ~ 120 L Q (K kms -1 pc 2 ) -1 
is within the range found in local ULIRGs (Solomon et al. 
1997) but lower than values estimated for dusty quasars 
(e.g. Carilli etal. l2fTo"2l Walter etal. HSDB) . All this im- 
plies that large amounts of molecular gas are less affected 
by the heating from the AGN and/or nuclear star for- 
mation than the region emitting in CO(7-6). A possible 
explanation for this finding is that the molecular disk ex- 
tends beyond the central region seen in CO (7-6) to a zone 
of more moderately excited molecular gas that dominates 
the global emission in C i and the low J CO lines. 
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